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I.  INTRODUCTORY  STATEMENT 

Previous  to  the  seventeenth  century  no  practicable  methods 
for  the  solution  of  equations  existed  which  are  -worthy  of  note.  Even 
in  1600  when  Vieta  extended  and  generalized  a  method  invented  by  early 
Hindu  mathematicians,  the  process  of  solution  was  so  laborious  that  a 
contemporary  mathematician  characterized  it  as  "work  unfit  for  a 
Christian."  *    Thus,  practically  all  the  important  work  in  the  solution 
of  equations  is  the  product  of  the  past  three  centuries. 

It  is  the  purpose  of  this  thesis  to  set  forth  the  important 
usuable  methods  for  the  solution  of  algebraic  and  transcendental 
equations.    Four  types  of  solution  will  be  considered:     literal,  numer- 
ical, graphical,  and  nomographic  methods.     It  might  be  well  to  mention 
that  in  addition  to  these  four  methods,  equations  may  also  be  solved 
by  mechanical  devices.    Many  ingenious  machines  have  been  invented  but 
their  cost  and  complexity  prohibit  their  extensive  use. 

II.  LITERAL  METHODS  OF  SOLUTION 

Literal  methods  of  solution  consist  of  formulae  for  the 
roots  of  an  equation  in  terms  of  its  coefficients.    By  this  method 
solutions  of  the  general  cubic  and  quartic  are  obtained.     It  has  been 
proved  by  Abel  that  there  exists  in  general  no  such  solution  for 
equations  of  higher  degree  than  the  quartic.** 

*    The  Calculus  of  Observations  -  ".Vhittaker  and  Robinson,  Page  ?9 
**  Source  Book  in  Mathematics  -  D.  E.  Smith,  Page  261 
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1.    Solution  of  the  Cubic. 

(a)     Tartaglia's  Solution  * 

Tartaglia's  solution  of  the  general  cubic  equation, 

x3  -f  b::*+  cx  f  d  =  0.  (l) 

In  (1)  substitute  x=-y  -  —  (2) 

o 

We  then  get  y3  j-  py  +  qsO.  .  (3) 

This  contains  no  term  in  yl  and  is  called  the  reduced  cubic. 

In  (3)  substitute  y  —  z  -  —  (4) 

3z 

Then  z3  f  q=0 

or  z    f-  qz    -  i—  =0.  (5) 

27  ' 

This  is  a  quadratic  in  z  and  oan  be  solved  by  mears  of  the  quadratic 
formula. 

If  we  let  R  s  5l  +  JL 

4    T  27 

Then  z3  =  -  |  +.y?7  (6) 

And  z'  =   -  I  -p".  (7) 

In  equation  (6)  there  are  three  cube  roots  zf  ,  z^,  and  z  .    We  then 
substitute  these  values  in  (4)  and  get 

Again,  substituting  in  (2)  we  find  the  roots  of  the  general  cubic 
(1)  to  be 

x   -  z    -  P  _  L  •    x-z-P  L  .    *  -  7  _  P    _  L                     r  frt 

*'       *'     3%"'     X*~          3^     3'    XJ-Z3  3l3-  3  * 

It  can  be  proved  that  the  three  cube  roots  of  (7),  z„,  z  ,  z  are 

T         3  x> 

identical  with  the  results  obtained  by  solving  (6),  hence  the  solutions 
given  by  (8)  constitute  a  complete  solution  of  (l). 


*    College  Algebra  -  W.  L.  Hart,  Page  248. 
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(b)     Trigor.ometr ic  Solution  of  the  Irreducible  Case.  * 

When  the  roots  of  a  cubic  equation  are  all  real  and  distinct, 
R  (see  preceding  section)  is  negative  so  that  the  formulae  used  in 
Tartaglia's  solution  yield  values  in  a  form  involving  cube  roots  of 
imaginaries.    This  is  called  the  irreducible  case.**For  this  we 
obtain  a  direct  solution  based  upon  the  well-known  trigonometric  identity 

cos3xs4cos  x  -  3cosx 
This  may  be  written  in  the  form 

z3  -  %Z  -  ^-cos3x=0.  (z=cosx). 
To  transform  the  reduced  cubic  into  this  one,  set  ysnz.  Thus 

The  two  cubic  equations  are  identical  if 


\  3      »      cos3x  =■  — £   —  \j-^3 


Since  R<0,  p<0  and  the  value  of  cos3x  is  real  and  numerically  <C  1. 
Hence  we  can  find  3x  from  a  table  of  cosines.  The  three  values  of  z 
are  then 

cosx,  cos  (x  t  120°),  cos  (x  t  24(5). 
Lultip  lying  these  by  n,  we  get  the  three  roots  of  the  reduced  cubic. 

There  exists  another  solution  of  the  cubic  known  as  the 
trisection  method,  which  is  applicable  whatever  the  nature  of  the  roots 
of  an  equation,  and  yrhich  uses  the  hyperbolic  functions  instead  of  the 
tr igonometr ic. *** 


*    Elementary  Theory  of  Equations  -    L.  E.  Dickson,  Pages  34,  36,  37 

**  This  can  be  done,  however,  by  Tartaglia's  Lethod. 

***The  Calculus  of  Observations  -  '"hittaker  and  Robinson,  Pages  124-125. 


c 


♦ 


2.    Solution  of  the  Quartic  (Ferrari's  Solution).  * 
Ferrari's  solution  of  the  general  quartic 

x^  \  bx3  +  ex3"  +  dx  +  esO,  (l) 
From  (1)  xH  +  bx3  =.  -  cx~  -  dx  -  e.  (2) 

If  we  add  ^bXx  tc  both  sides  of  (2)  the  left  hand  member  becomes  a 
perfect  square,  and  we  have 

(xx  +  jjtac)*  =  (ib*  -  c)x*-dx  -  e.  (3) 
Let  y  be  a  number  to  be  determined  later.     On  adding  |V(x*"  +  ^bx)  +  4y*"J 
to  both  sides  of  (5)  we  obtain 

(xl  +  Jbx  +  §yf  =  (Jb*  -  c  +  y)x*  +  (Jby  -  d)  x  +(j  y*  -  e) .  (4) 
The  right  side  of  (4)  is  a  quadratic  in  x  which  will  be  a  perfect  square 
if  the  discriminant  is  zero,  i0e.,  if 

(Jby  -  df  -  4(^b*-  -  c  ^(hr*-  e)  =  0.  (5) 
Yftiich  simplified  becomes 

y3  -  cy*  +  (bd  -  4e)y  -  bze  +  4ce  -  d*=  0.  (6) 
This  is  called  the  resolvent  cubic  for  equation  (l)  and  can  be  solved 
by  either  Tartaglia's  solution  or  the  trisection  method. 
Let  y  «  s  be  any  root  of  (S),  and  "substitute  in  (4). 
The  right  side  of  (4)  becomes  the  square  of  some  linear  function, 


(hx  +  k) ,  and  (4)  becomes 

(x*  +  Jbx  +  Js)*  =  (hx  +  fcf  (7) 

Then  xx  +  §bx  +  ths  -  hx  +  k  (8) 

or  x2-  +  Jbx  +  js  =•  -(hx  +  k)  (9) 


Both  (8)  and  (9)  can  now  be  solved  by  use  of  the  quadratic  formula.  The 
four  values  of  x  thus  obtained  are  the  roots  of  (l)  . 

It  is  interesting  to  note  that  by  the  use  of  the  formulae 

*    College  Algebra  -  W.  L.  Hart,  Page  249 
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derived  for  ""artagli?.  's  solution  of  the  general  cubic  we  could  write  • 
a  formula  for  the  root  g  of  the  resolvent  cubic  in  terns  of  the 
coefficients  of  (l).    We  could  then  obtain  formulae  for  h  and  k  and 
hence  for  the  four  roots  of  the  general  quartic  in  terms  of  its 
coefficients.     These  formulae,  however,  would  be  too  complicated  for 
practical  use. 

3.    Series  Formula  for  the  Root.  * 

This  method  furnishes  a  literal  formula  in  which  the  root 
is  obtained  by  mere  substitution.** 

The  root  of  the  equation 


X  1 

3 


0»  a#  *  af  *     <  {  +  a„xY  f   


v;hich  is  the  smallest  in  absolute  value,  is  givsn  by  the  series 


•  o  i. 
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For  a  proof  of  this  formula,  let  the  roots  of  equation  (l) 
be  x  ,  x a  x_  ,   ,  X^. .     Then  if  z  is  any  number  whose  modulus  is 


smaller  than  each  of  the  moduli  of  the  roots,  we  have 


a     t  a,  z  *  a  z' 


 _  —  ...  ...  ~i     «■  » 


=  (1  -  z/x,  )(1  -  ZJ6J.....U  -  z/xn) 

i  -  z/x,     1  -  z/XjL  1  -  Z/X^ 

*    The  Calculus  of  Observations  -  Vtti  it  taker  and  Robinson,  Page  120 
**    It  will  be  noticed  that  this  method  is  applicable  to  equation*  of 
the  nth  degree  in  contrast  to  the  other  literal  methods. 
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< 


=  (1  t  z/x,  f  zZ/x*f +  zAl+  z  /\  +  )  +  z/xn+  z  /-V--..) 


=      +    Pf  2  ^  S*    +  P3  Z  * 


where  P^  denotes  the  sum  of  the  homogeneous  powers  and  products  of  the 
reciprocals  of  the  roots  taken  r  at  a  time. 

Therefore      a0=  (a6  h  a|Z  +  a^z*   ^   )(l  f  P#.a  t  P^z^  f  P3  z3  +....) 

Equating  coefficients  of  powers  of  z  we  have 
0  .  a,   t  a.?,, 
C  =  ax  f      r(    t  a0Px  , 
=     ,  / 


P    +  a 

f 

a 

0 

p 

3 

• 
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1 
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0 

0 
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a3 

Z 

0 

'© 

o 

/ 

3 

itc. 

Now  since    r  *1sa2  -  a^a.  .  we  see  that  tho  first  two  terms  of  the 
series  (2)  are  equivalent  to  the  single  term 
  atf  ay  P| 

I  •  ,1 


(3) 


or 


(4) 


Similarly  vre  can  show  that  the  sum  of  the  first  three  terms  of  the  series 

P  P 
is  equal  to  <p*  ;  of  the  first  four  terms,  wL,  eto«     Thus  we  see  that  the 

\3  *V 


sum  of  the  first  s  terms  of  the  series  is  equal  to 


5-,  . 


S 


Let  us  for  simplicity  consid3r  the  case  when  n  =  2 ,  where  x    has  a  smaller 
modulus  than  the  root  x^ .     Thou,  recalling  the  definition  previously 
have 


given  of  P^  ,  we 


i-i 


--  S-T. 

-  I  -"-x 


_  3-3 


5-/ 


4_  +  + 


i-   + 
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Multiplying  numerator  ard  denominator  by  x    and  factoring  out  xf  ,  we  get  , 

x,       xf  X?" 

S-l  - 


-    X,  .                    --t                      5-/  ^5 
1   +-  — —  -4-  -^i—  +-       .    X  ^   4-  it. 


x    —  +  rr  t-  ...4-  ^r-  +  x, 

and  since  the  numerator  and  denominator  are  the  same  with  the  exception  of 
s 

the  term  in  the  denominator,  and  since    ~|<  1,  this  gives  at  once 

T  ^r        S  -l 

Similar  reasoning  leads  to  the  result  when  n>  2. 

The  series  converges  rapidly  when  the  ratio  of  the  smallest 
root  to  erery  one  of  the  other  roots  is  small.     In  calculating  any  root 
of  a  given  equation  by  the  formula  it  is  often  advisable  to  transform 
the  equation  by  two  or  three  root  squaring  operations;*  or  if  an 
approximate  value  of  the  root  is  known,  to  reduce  the  roots  of  the 
given  equation  by  a,  where  a  is  an  approximate  value  of  the  required 
root.     In  either  event  the  required  root  of  the  new  equation  will  be 
small  as  compared  with  any  of  the  other  roots. 

As  a  numerical  example  consider  the  equation 
x3  -  2x  -  5  =  0, 

whore  we  are  to  find  the  root  which  lies  between  2  and  3.    Since  2 
is  an  approximation  of  the  root  we  will  first  transform  the  given 
equation  into  one  whose  roots  are  decreased  by  -2.     The  transformed  equation 
becomes  x    f  6x    f  lOx  -  1=  0.  ** 

1  a  a0s  - 1 ,  a  =  10 ,  a4  =  6 ,  1,3  1, 

/    z  31 


/    Z    3  1 

|a  l^Sl°6'  a«a'a*k121> 


*  See  root-squaring  process,  page  19. 
**  See  Ruff iai-Kornsr  Method,  page  9. 
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X    3  V 

'  A  J 
«     /  A 


*    x  J  v 

0     1    V  J 

.    /  x 


=  11,356. 


The  smallest  root  of  the  equation  is  therefore 


26 


95 


10      1O10S  '  105«1121  1121-11856 
.1  -  .005660  +  .0002188  -  .0000072  a  .0945515 
Thus  the  root  of  the  given  equation  is  2.0945515. 


III.     NUMERICAL  METHODS  OF  SOLUTION. 

The  process  of  the  numerical  method  of  solution  consists 
mainly  of  arithmetical  computation.     It  is  on  the  whole  the  most 
useful  type  of  solution,  especially  when  a  high  degree  of  accuracy  is 
des  ir  ed. 

1.     Equations  in  One  Unknov<-n. 

(a)     The  Ruffini- Horner  Ifothod. 

The  Ruf f ini-Korner  method  of  finding  the  roots  of  an  equation 
is  a  less  laborious  extension  of  the  method  of  Vieta,  and  should  be 
familiar  to  all  who  have  taken  a  course  in  college  algebra.  After 
locating  the  roots  between  consecutive  integers,  an  approximation  may 
be  formed  to  almost  any  desired  accuracy  be  decreasing  the  roots  of 
the  equation  by  means  of  successive  transformations.    Given  the  equation 

x3  -  2x  -  5=  0. 

It  may  be  shovm  by  graphirg  that  one  root  of  the  equation  lies  between 
2  and  3,  so  we  first  decrease  the  roots  of  the  given  equation  by  2. 


-  3  - 


6 


- 


2 

_2_ 
4 
2 


2 
8 
10 


-  1 


This  gives  the  transformed  equation  x3  +  6x    *  lOx  -  1=0. 
Equating  the  last  two  terms  of  this  equation  to  zero  for  an  approximation 
we  find  x  =  .l  to  be  an  approximate  root.    By  substitution  or  synthetic 
division  there  is  found  to  be  a  change  in  signs  of  f(x)  between  x  =  .09 
and  x  =  .10,  (i.e.  the  curve  cuts  the  x-axis  between  x  =  .09  and  x=  .10). 
Thus  we  next  decrease  the  equation  by  .09. 

.09 


1 

6 

10 

-  1 

0 

.09 

.548 

.94932 

1 

6.09 

10.548 

-  .05068 

0 

.09 

.556 

1 

6.18 

11.104 

0 

.09 

1 

6.27 

is 

repeated  until  the  desired  aocuracy 

1 

6.27 

11.104 

-  0.05068 

0 

.004 

.025 
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1 

6.274 

11.129 

.00616 

0 

.004 

.025 

1 

6.278 

11.154 

0 

.004 

1 

6.282 

1 

6.282 

11.154 

-  .00615 
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.0005 

.0031 

.00558 

1 

6.2825 

11.1571 

.00058 
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.0005 

.0031 

I  .004 


.0005 


1 

0 


6  .  2830 
.0005 


11.1602 
1  6.2835 

Equating  the  last  two  terms  of  this  transformed  equation  to  zero  we  find 
x  =.00005  to  be  the  next  approximation. 
Thus  the  root  is  2.09455. 

(b)    The  Method  of  Iteration.* 

When  a  numerical  equation  f (x)  =-  0  can  be  expressed  in  the 


form 
* 


U) 


The  Calculus  of  Observations  -  Whittak^r  and  Robinson,  Pages  81  -  83. 
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the  real  roots  may  he  obtained  "by  the  process  of  iteration.    We  find 
by  means  of  a  graph  or  some  other  method  an  approximation,  xo ,  of  the 
desired  root.     This  value  is  then  substituted  in  the  right  hand  member 
of  (1)  and  a  new  approximation  x(    is  given  by  the  equation 

x,  =  f(*«)« 

Substituting  this  new  value  x.  we  get 

xi=  ^(x,  ). 

Succeeding  approximations  are  found  in  a  like  manner, 

x?=  Mx  ) 
x¥=  ^(x3) 

x„  =  ). 

An  interesting  fact  in  connection  with  the  method  of  iteration 
is  that  a  mistake  in  the  numorical  calculation  does  not  invalidate  the 
entire  result.     If,  for  example,  a  mistake  were  made  in  calculating  x.t 
from  xft,  the  erroneous  value  %  might  have  been  obtained  correctly  by 
starting  from  a  different  first  approximation  x^ ;  and  since  xe  is 
arbitrarily  chosen,  the  true  solution  ma;/  be  reached  by  way  of  X,  as 
well  as  by  way  of  r.f  . 

The  method  of  iteration  may  easily  be  illustrated 
geometrically.     If  the  function  f(x)— 0  is  written  in  the  form 
f/(x)  =  fJt(x)  and  the  respective  curves  are  drawn,  the  abscissae  of  the 
points  of  intersection  of  these  two  curves  are  evidently  the  real  roots 
of  the  function. 

An  iterative  process  for  finding  them  may  be  obtained  as 
follows:     select  any  point  xQ  on  the  x-axis  in  the  neighborhood  of 
the  abscissa  of  the  point  of  intersection  of  the  curves.    From  x 

O 


-10- 
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draw  a  straight  line  parallel  to  the  y-axis  until  it  cuts  the  curve 
v.-hich  has  the  slope  of  lesser  magnitude .    Then  draw  a  line  parallel 
to  the  x-axis  until  it  meets  the  second  curve.    From  this  point  draw 
another  parallel  to  the  y-axis  until  it  again  meets  the  first  curve, 
and  continue  this  process. 


A.  B. 


The  accompanying  diagrams  *  illustrate  the  two  main  types 

resulting  from  this  process,  according  as  the  slopes  of  the  two  curves 

have  the  same  or  different  signs  for  the  abscissae  x  ,  x  .  x   

In  diagram  £  the  abscissae  x  ,  x  ,  x  ,  .......  are  all  on  the  same  side 

°       '  2. 

of  the  root  x  and  the  lines  approach  the  point  of  intersection  in  the 
form  of  a  staircase.     In  diagram  B  the  lines  approach  the  intersections 
spirally.     In  the  first  case  the  derivatives  of  the  curves  have  the 
same  sign  near  the  point  of  intersection  and  in  the  second  case  the 
signs  are  unlike. 


*    The  Calculus  of  Observations  -  7,'hittaker  and  Robinson,  Page  32. 
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Consider,  again,  the  equation 

x3  -  2x  -  5-0. 
We  first  -write  this  in  the  form. 


'1 


2x  f  5 


Then  if  x,  =•  2  we  find 


=  2.025, 


and 


•\  2  2.025  +  5  _ 
Vl|     2i 02 5      "  -.j-^j-. 


2.0802, 

4"  = 

2.0979, 

6  = 

2.0936, 

X7  = 

2.0948, 

2.09457. 

(c)     The  Method  of  False  Position  (iiegula  Falsi)  * 

The  oldest  method  for  computing  the  real  roots  of  a  numerical 
equation  is  the  method  of  false  position  or  "regula  falsi".     This  is  the 


*    Numerical  Mathematical  Analysis  -  J.  B.  Scarborough,  Pages  174-176. 
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numerical  equivalent  of  the  graphical  method  of  repeated  plotting.  We 
find  two  numbers  xf  and  xz  between  which  the  root  lies.    Since  the  curve 
must  cross  the  axis  between  these  two  values  then  y7  and  yfc  must  have 
opposite  signs.     A  smooth  carve  will  be  practically  sbraight  for  a  short 
distance,  hence  it  maybe  assumed  that  the  change  in  f(x)  is  proportional 
to  the  change  in  x  over  a  short  interval,  as  in  the  case  of  linear 
interpolation  from  logar ithmic  and  trigonometric  tables.     The  method 
of  false  position  is  based  on  this  principal,  for  it  assumes  that  the 
function  f(x)-=0,  is  a  straight  line  between  the  points  (x,  ,  y,  )  and 
(xz,  yz ) ,  these  points  being  on  opposite  sides  of  the  x-axis. 

•  The  accompanying  diagram  represents  a  magnified  view  of 
that  portion  of  the  graph  of  f(x)  —  0  between  (x, ,  y, )  and  (x»,  ya.). 


P(*.;y-) 


Then  from  similar  tr  i  angles  FIG  and  PRQ    we  hav< 

h         xx  -  x, 


The  value  of  the  desired  root  then  is 

x  =  x,   +       —        *■  h. 
Hence  x  =  X,  +  3  M  . 

\y,\  * 
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This  value  of  x  is  not,  hoover,  the  true  value  of  the  root 
because  the  graph  of  the  function  f(x)  =  0  is  not  a  perfectly  straight 
line  between  P  and  Q.     It  is  merely  an  approximation  of  the  true  root. 
In  the  application  of  tho  regula  falsi  method  we  compute  a  short  table 
of  corresponding  values  of  x  and  f(x)  for  equidistant  values  of  x 
(units,  tenths,  hundredths,  etc.)     Then    by  the  use -of  the  formula 

hi  "W 

v/e  compute  corrections  to  be  applied  to  the  previously  obtained 
approximations.  To  make  this  method  clear  the  following  example  is 
introduced. 

Example;     Compute  the  real  root  of  x  Iog/Ox  -  1.2=  0.* 
Having  previously  determined  that  the  root  lies  between  2  and  3,  we 
proceed  with  the  following  table,  using  logarithms  to  more  decimal 
places  with  each  succeeding  approximation, 
x  y 


1st 

2 

-0.6 

- 

app. 

3 

+0.23 

tiff. 

1 

0.83 

x,  = 

2nd 

2.7 

-0.04 

hJL  = 

app. 

2.8 

+0.05 

diff . 

.1 

0.03 

TC  — 

3rd 

2.74 

-0.0006 

h3  = 

app. 

2.75 

+0.0081 

diff. 

0.01 

0.0087 

x3  = 

4  th 

2.7406 

-0.000059 

h*  = 

app. 

2.7407 

+0.000045 

diff. 

0.0001 

0.000084 

1  •  0.6 


0.83 
2.72 

0.1  •  0.04 


-  0.72 


0.09 
2.74 

0.01  •  0.0006 

0.0087 
2.7407 


z.  0.044 


= 0.0007 


0.0001-0.000039 

0.000084 
2.74065 


-0.000046 


*    Taken  from 

numerical  Mathematical  Analysis,  J.  B.  Scarborough,  Page  175 
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e 


c 


Another  method  of  deriving  the  rule  of  false  position  is 
as  follows.  Let  f  (x)  "be  any  curve  cutting  the  x-axis  and  (x,  ,  y,  ) 
and  (x^ ,  yz )  he  two  points  on  this  curve. 


The  equation  of  the  chord  correcting  these  points  is 

y  -  y*  -  Vjl  -  :'«  , 

X  -  x,       xx  -  X, 

If  x  is  to  he  a  root  then  y=0,  hence  we  have 

-  y>    _  y*  -  y# 

X    "    X/  XJL    "  X/ 

«  „  ■       U*.  -  x/)y> 

and  x=-(.  . 

'       y*  -  y> 

By  use  of  this  formula  let  us  solve  the  equation 

x3  -  2x  -  5~  0. 

Here  we  can  take    x'-  2    _  (5  ~  2)(-  1)  -  2  f  L_  =  2.06. 

16+1  17 


How  let  us  take  x7  =  2.06,  x&=  2. 10, 


;hen 


2#06  _  (2.10  -  2.06)  (-  .5  7  8184) 


.062  t  .378184 

=  2.06  f  .034=  2.094 
Again  if  we  let  x  =2.094,  and  x  =  2.096,  since  f  (2.096)  or  ya  is 
po  s  it  iv  e ,  -we  have 

XW_    2  0o4       (2.096  -  2.094)  (-.006155416) 
X    S      "  .016180736  f  .006153416 

=  2.094  f  .0005510  =  2.0945510 
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(d)    The  Newton-Raphson  Method. 

When  the  derivative  of  f(x)  is  a  simple  expression  and  easily 
found,  the  real  roots  of  f  (x)  =  0  can  he  computed  rapidly  "by  a  process 
called  the  Newton-Raphson  method.     The  underlying  idea  of  the  method 
is  due  to  Newton  who  calculated  a  set  of  successive  equations  v/hose 
roots  were  the  successive  remainders  between  the  approximation  and  the 
true  value  of  the  root,* 

For  example,  to  find  the  root  "between  2  and  3  of 

x3  -  2x  -  5  =  0 
Newton  replaced  x  bv  2  t  P  and  obtained 

p3  f  Sp*  f  lOp  -  1*  0. 

Since  p  is  a  decimal,  he  neglected  the  first  two  terms  and  set  lOp  -  1= 

so  that  p=0.1  approximately.    Replacing  p  by  0.1  f  q  in  the  preceding 

cubic  equation,  he  obtained 

q3   +  6.3qZ  f  11.23q  +  0.061=0. 

Dividing  -0.061  by  11.23,  he  obtained  -0.0054  as  the  approximate  value 

3 

of  q.     ileglecting  q    and  replacing  q  by  -0.0054  +  r,  he  obtained 

6.3r2  f  11.16196r  t  0.000541708=0 
Dropping  6.3r*,  he  found  r  and  hence 

x=  2  f  0.1  -  0.0054  -  0.00004853 
=2.09455147.  ** 
However,  the  method  as  now  used  is  due  to  Raphson.  Instead 
of  iterating  towards  a  root  in  rectangular  steps  as  in  the  method 
of  iteration,  -we  now  proceed  in  oblique  steps.     Using  rectangular 
coordinates  consider  the  graph  of  y  =  f(x)  and  the  point  P  on  it  with 

*    Numerical  Mathematical  Analysis  -  J.B.  Scarborough,  Page  178 
**  Flementary  Theory  of  Equations  -  L.  '.  Dixon,  Page  109. 
Note  similarity  to  Ruff ini-Horner  Method,  page  9. 
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the  abscissa  OQsa 


q      t  r, 

Let  the  tangent  at  P  meet  the  x-axis  at  T  and  let  the  graph  meet 

the  x-axis  at  S.     Take  h  =  QT,  then 

QP=f  (a),  f '  ( a)  =  tanXTP^  _       =  _  £iil  . 

.       QT  h 
.  h  _     f (a) 

In  the  graph,  OT  =  a  +  h  is  a  better  approx imat ion  of  the  root  OS 
than  the  value  OQ  —  a.     To  get  a  still  more  accurate  approximation 
repeat  this  process  (as  indicated  in  the  diagram  "by  the  dotted  line). 
This  may  he  continued  until  the  desired  accuracy  is  obtained. 

We  may  also  derive  the  same  formula  by  the  use  of  Taylor's  theorem. 
If  a  is  an  approximate  value  of  the  root  of  f(x)=0  and  h  is  the 
correction  which  must  be  applied  to  give  the  exact  value  of  the  root,  then 

x  =  a  f  h. 
The  equation  f  (x)  =0  then  becomes 

f(a  f  h)  =  0. 
Expanding  this  by  Taylor's  formula  we  have 

f(a  +  h)  =  f(a)  +  f'(a)h  +  f''(a)|\-  


Since  h  is  a  very  small  quantity  we  disregard  the  higher  powers  of  h 
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and  take  only 

f(a)  f  f'(a)h=0, 

f(  a1) 

The  improved  value  of  the  root  is  then 

f(a,) 

Similarly,  »,  *  h*~  a/~  f'U,  ) 


"ev/ton's  formula  also  comes  directly  from  the  formula  for  the  "Regula 
Falsi"  method.  If 


A 

x  =  x ,  —  —  i 


and  x  =  x 


/  a 


/\v  dv 

Then  as^x-»  0,  "becomes  (i.e.  -we  replace  the  chord  v;ith  a 

tangent)  and  we  have 

x-  x    —  **- 


Solving  again  for  the  root  of  the  equation 

xJ  -  2x  -  5=0 


vrhich  lies  between  2  and  3  we  have 


■ 


°v     f7!^)"  2       10    -2  +  1Q-  2.1 

CGI 

j=  i- •  1      11  ^  3      2.1  ~~  »00o4^  2.0946. 
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2.0946  -  t222||^g  -  2.0946  -  .000047=  2.094563 
a.=  2.094563  -  TT^TgTTO  =  2.094565  -  .000007  =  2.094556 

5  XX  .  ID  J-OOtz 

(e)    Graeffe's  Root-Squar ing  Method. 

The  numerical  methods  of  solution  thus  far  considered  are 

♦ 

applicable  only  in  solving  for  real  roots  of  equations.     The  root- 
squaring  method  of  Graeffe  can  be  applied  to  either  real  or  complex 
roots . * 

The  underlying  principle  is  this:     the  given  equation  is 
transformed  into  one  whose  roots  are  high  powers  of  the  roots  of  the 
original  equation.     The  roots  of  this  transformed  equation  are  widely 
separated  and  because  of  this  fact  are  easily  found.    For  example,  if 
two  of  the  roots  of  a  given  equation  are  3  and  2,  the  corresponding 

^77  777 

roots  of  the  transformed  equation  are  3    and  2     ,  where  m  is  the 

power  to  which  the  roots  of  the  given  equation  have  been  raised. 

Thus,  if  m=  64,  we  have  3    -    10  ,  2    =  10  .     The  Wo  roots 

of  the  given  equation  were  of  the  same  order  of  magnitude  but  in  the 

transformed  equation  the  larger  root  is  more  than  a  hundred  billion 

times  as  large  as  the  smaller  one.     The  ratio  of  the  roots  of  the 

2 

given  equation  is    "%  »  but  in  the  transformed  equation  it  is 


Dandelin  and  Lobachevsky  also  developed  this  method  independently, 
antedating  Graeffe. 
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-^—so.iiir~  ~  iQ  /t.tf-  '  or  JZ?  <  0.000000( 

tra.nsforr.3d  aquation  is  therefore  negligible  in  comparison  with  the 
larger  one.     The  roots  of  the  transformed  equation  are  said  to  he 
separated  when  the  ratio  of  any  root  to  the  next  larger  is  negligible 
in  comparison  with  unity. 

The  transformed  equation  is  obtained  hy  repeated  applications 
of  a  root-squaring  process.     The  first  equation  is  transformed  into 
one  whose  roots  are  the  squares  of  the  roots  of  the  given  equation. 
Then  this  process  is  repeated  with  the  second  equation  and  a  third 
equation  is  obtained  in  which  the  roots  are  the  squares  of  the  roots 
of  the  second  equation,  and  hence  the  fourth  powers  of  the  roots  of 
the  original  equation.     This  process  is  continued  until  the  roots  of 
the  last  transformed  equation  are  completely  separated. 
Let  the  given  equation  be 

f(x)  =  yx"  ♦  a,x"-'   f  Vc"-*   ,    *  *  V  0.  (1) 

Then  If  x  ,  x  ,   x    be  the  roots  of  this  equation  we  can  write  it  in 

the  equivalent  form 

f(x)=  a0(x  -  x,)(x  -  X<f)(x  -  x^)  (x  -  0.  (2) 

How  let  us  multiply  (2)  by  the  function 

(-if  f(-x)  =  (-l)%0(-x  -x,)(-x  -  xj  (-x  -  (3) 

c  ao(x  t  x(  )  (x  f  xz)  (x  +  x^  ) . 

The  result  is 

(-l)rtf(-x)f(x)=a*(xt  -  x*)(x*  -  x*)  (x1  -  x*).  (4) 

Let  x  =  y.     Then  (4)  becomes 

/(v)  =  %(V  -  xf)(y  -   Cy  "       =  0.  (5) 
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The  roots  of  this  equation  are  x*,  x.t  ....  x*and  are  thus  the  squares 
of  the  roots  of  the  given  equation  (l).    Hence  to  form  an  equation  whose 
roots  are  the  squares  of  those  of  f(x)  =  0,  we  merely  multiply  f  (x)=  0 
by  (-lff(-x).  * 

This  multiplication  can  be  carried  out  in  a  simple  routine 
manner.     Consider  the  sixth  degree  equation. 

f(x)  -  aox'  4-  s,  yf  f  a^x*  f  ^  x3   f  a^x**  f  afx  *  afc=0. 
Then        (-if  f  (-x)  —  aex6  -  a,  rf  f-  syx*  -  a3x3   i  a^x*"-  ajc  i»-  a^. 
By  actual  multiplication  we  find 


(-l)6f(-x)f(x)=  aV*  -        aV°  t 


S-  x  6 

+  2a0a4|       -  2a,  aJ      t  2aJLa), 
f  Ea^j      -  2a,  a, 


f  2a<}at 


"  2a,  a, 
f  2aAat 


<■  2av  2^ 


f  0.  (6) 


The  law  of  formation  of  coefficients  in  (6)  (and  likewise  for 
eqaations  of  any  degree)  may  therefore  be  stated  thus:    the  coefficient 
of  any  power  of  x  is  formed  by  adding  to  the  square  of  the  corresponding 
coefficient  in  the  original  equation  the  doubled  product  of  every  pair 
of  coefficients  which  stand  equally  far  from  it  on  either  side,  these 
products  being  taken  with  signs  alternately  negative  and  positive.  ** 
Case  I.    Roots  all  Real  and  Unequal. 

Since  the  relations  between  the  roots  xf  ,  x%,  ....  x^and 

^  )  coefficients  ao,  a.f  ,   \°-  Ssnera-1  equation  of  the  nth  degree 

a0x    t  a,  x        f    *  o.fj^x  t  a^=.  0 

are  =s.-i  fet  +  x.  +    +  x  ), 


0 


*  Numerical  Mathematical  Analysis  -  J.  b.  Scarborough,  Page  199 
**  The  Calculus  of  Observations  -  Whittaker  &  rlobinson,  Page  103 
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< 


+  (x,  XjL  +  x(  x  t 


o 


(-1)  x,*^  a 


it  follows  that  the  roots  xf,  x*",   x^and  the  coefficients  b_  ,  b.  , 

 b  of  the  final  transformed  equation 

b.  (x-f    *  b,  (ff~'+   *  b„_f  x"  +  b„  =  0 

are  connected  by  the  corresponding  relations 


-  x7(i  +      +   +ia  ) 

7*i         **T  ^7 


X7"  X1 

/  / 


b*.  _    *   m       m  .m  »»   m,         x3       x-  ,  Xj  , 


x/ 

b»  /  m  «i  m  t»i  -jn  ^  m   m  mt  ™ 

-jL.  -    -    ( X  X.  X      +  X  X  X      f  )  -  -  X  X  X    (  1  4-         4.   .  .  .  .  .i.  ±JL  4. 

b,  »X3r      /    XV  ;  /    a.    3  v  X  m       •  •  •  '  T  £^  +-    .  .  ,  ) 


hi 

bo 


Heir  if  the  order  of  magnitude  of  the  roots  is 

N>l*J>hl  •••••>tei  » 

it  is  evident  that  -when  the  roots  are  sufficiently  separated  the  ratios 


x^      x™ .  etc.  are  negligible  in  comparison  with  unity,    Kence  the 

r»i    '    -  tn  ' 

relations  between  roots  and  coefficients  in  the  final  transformed 


e  qua.  t  ions  are 


b/  ">    b~  m    b»  m  m  -m 

-  -r        •     _  -r  v         3    _     -x  x  3 


bo  '  *  3  * 

Dividing  each  of  these  equations  after  the  first  by  the  preceding  equation, 

we  obtain 


b7"  ~  -k*  bT-    -3'  sj_ 
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< 


'lence  from  these  and  the  equation  xf  we  get  * 

The  root-squaring  process  has  thus  broken  up  the  original 
equation  into  n  simple  equations  from  which  the  desired  roots  can  be 
found  with  ease. 

The  question  arises  as  to  how  many  root-squar ings  are  necessary 
to  break  up  the  original  equation  into  linear  fragments.     This  cannot  be 
determined  beforehand,  and  in  practice  the  root-sq>3ar ing  method  is 
continued  until  the  doubled  products  in  the  second  row  have  no  effect  on 
the  coefficients  of  the  next  transformed  equation. 

Since  the  coefficients  in  the  given  equation  are  not  in  general 
all  positive,  the  signs  of  the  doubled  products  will  not  occur  in 
regular  order  as  in  the  literal  equation  used  to  illustrate  the  process. 
Hence  care  must  be  taken  to  avoid  the  possibility  of  making  a  mistake 
in  s  igns . 

"xample:    Find  all  the  roots  of  the  equation 

Xr  -  3x    -  x    f  4x  f  2=0, 
The  preliminary  work  of  separating  is  shown  on  the  following  page. 

It  is  evident  that  further  squaring  will  simply  give  the 
squares  of  the  coefficients  in  the  last  line  of  the  table,   and  we 
trierefore  stop  with  the  64th  powers  of  the  roots.     Thus  we  have  the 
four  simple  equations: 

zy   -  3.617  •  10  =0. 

(-8.617  .  1027)x^   f  2.043  *  10W=  0. 
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I  III 


1    I  I 


log 

8.617- 

ioA7 

.43649. 

64 

10*  - 

io" 

2.045- 

log  6.317- 

Ic  ; 

4.354  - 

J* 

10 

log  2.043- 

¥/ 

10 

1.845- 

10"  - 

log  4.354- 

3X 

10 

•   H'        -  4.:    •   3*=  o. 

(-4.354-10    )x*7  f  1.845-10  =  0. 
Solving  these  "by  logarithms,  we  hav 
log  x,= 

log  x2=    ~"   —  — ~-  — — ~ — -  .6' 

log  4.354-10*    -  log  2.043- 10* 
log  x3= —  ~   9.86451  -  10. 

log  1.845-10"    -  log  4.354-10*" 
log   ~  ~  9.79105  -  10. 

Thus  we  have        x#s*  2.732,  x^s  1.618,  .732,  x^s  .618. 

The  signs  of  these  roots  are  yet  to  he  determined.     To  do  this  we  first 

apply  Descartes'  rule  of  signs  and  find  that  there  cannot  he  more  than 

two  positive  roots,  nor  more  than  two  negative  roots.     Then  we  substitute 

in  the  given  equation  "both  positive  and  negative  approximations  of  the 

roots  to  see  which  value  comes  nearer  to  satisfying  the  equation.  In 

mt  the  r. 

x i  — -  2 .  id 2 

xz-  1.618 

v-  -T5Z 

xs-  .618 
V 

Case  II.     Complex  roots. 

Then  some  of  the  rooL-s  of  an  algehraic  equation  are  complex, 
the  equation  cannot  he  expressed  as  a  product  of  linear  factors  with 
real  coefficients;     Such  an  equation  can,  however,  he  expressed  as  the 
product  of  real  linear  and  real  quadratic  factors,  each  quadratic  factor 


c 


corres  pond  tag  to  a  pair  of  complex  roots.     The  root-squaring  process 
can  therefore  never  break  up  such  an  equation  into  linear  fragments 
as  is  the  case  vfoen  all  the  roots  are  real  and  distinct.    When  an  equation 
has  complex  roots,  the  root-squaring  process  always  breaks  it  up  into 
linear  and  quadratic  fragments.     The  real  roots,  if  any,  are  found  from 
the  linear  fragments  as  in  Case  I,  while  the  complex  roots  are  found 
from  the  quadratic  fragments. 

In  transforming  an  equation  by  the  root-squar ing  process  the 
presence  of  complex  roots  is  revealed  in  two  ways:     (l)     the  doubled 
products  do  not  all  disappear  from  the  first  row,  and  (2)  the  signs 
of  some  of  the  coefficients  fluctuate  as  the  transformations  continue.* 

After  the  given  equation  has  been  broken  up  into  linear  and 
quadratic  fragments  we  can  find  the  complex  roots  by  solving  the 
resulting  quadratic  equations  for  x'Tn  and  then  extracting  the  m**  root 
by  means  of  DeT'oivre's  theorem.     However,  this  method  leads  to 
ambiguities  of  sign  in  the  computed  roots  which  are  not  easily  removed. 

When  there  exists  only  one  pair  of  complex  roots  in  an 
equation,  the  following  method* *is  applicable. 

Consider  the  equation 

Tl  71-/  77-.X- 

a0x    f  a,x        +  V:   *  »;,=  0 

where  the  real  roots  x7  ,  x^ ,  x^   x^,  and  the  square  of  the  modulus 

of  the  complex  roots,  r2-,  have  already  been  computed^ ** In  order  to  find 
the  complex  roots  re     ,  re     ,  we  write  them  in  the  form  u  +  iv  and 
u  -  iv,  where  u=rcos9,  v  =  rsinQ,  so  that  re      f  re     s2u,     Then  from 

the  relationship  of  the  roots  and  coefficients  of  an  equation, 

*  This  is  due  to  the  presence  of  the  cosine  which  changes  signs  as  the 

amplitudes  are  continually  doubled. 
**  The  Calculus  o?  Observation  -  V  hittaker  and  Robinson,  Pages  113-116. 
***  The  value  of  r*  is  obtained  from  the  quadratic  fragment. 
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we  have  x,    fx,  fx.   +   fx.   f  2u=  ' 

i         *•       j  i  a 

o 

from  which  u  may  be  found.     Since  u  and  r*  are  known,  then  v  may  now 
be  found  from  the  relationship  r*=  u*"  +■  t  ,     The  two  complex  roots  are 
thus  determined. 

•  Again  suppose  we  have  the  equation 

a0x    f  a/  x        f  fc^x   .+  a,,.,   x  &  a;?=-  0 

which  has  two  pairs  of  complex  roots  r  ete  ,  r  e~'      and  r^e'9 ,  r^eT10  . 
Let  us  assume  that  the  real  roots  x/  ,  x^,  x^  ,  .....  x^,  and  the  squares 
of  the  moduli  r*"a  and  rA   have  already  been  computed.     If  u^s  r,  cos9, 
v  =  rf  sin9,  uA=  r^cos^,  v^=-  r^sin^,  the  complex  roots  are  denoted  by 
u/     iv,  ,  and  u  -±  iv^  ,  and  their  sum  is  2(uy  f  u^).    We  may  write  the 
given  equation  in  the  form 

(x  -  xy  )(x  -  xz  )(x  -  x3  )  (x*  -2ux  1-  r/")(xx- 2ujLx  +  £)-0. 

Thus  we  see  that  the  coefficient  of  x  in  this  equation  is  a  linear  function 
of  u,  an  d  u^.    Hence  we  have 

Solving  this  equation  simultaneously  with 

x,    t  x^  *  X.  f    f  x^  f  2u,   y  2u^=  -  £l 

(obtained  as  before)  we,  determine  the  value  of  u,  and  u^,  and  then  y,  and 
v^  by  substitution  in  tte  equation         r,a  -  u7  ,  v^"=:  r*  -  u^.     A  similar 
method  may  be  applied  to    solve  equations  with  more  than  two  pairs  of 
complex  roots. 

A  similar  application  of  this  method  is  obtained  by  sub- 
stituting x  =1  in  the  general  equation 

y 

a„x    f-  a,  x        f  a^x        4-   *°n-i  x  *  °* 
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The  result  after  clearing  fractions  is 

a^y    t  i„.,  y       f  a„-^y       *  .**,y* 

We  then  have  the  two  relationships 

X,     *  X*.   *  X3    * «^x77=  -  ~ 

and  L  +  L  +   ^I=_JLiiri. 

X/       XA-      X3  X"  a  ^ 

If  we  apply  this  to  the  preceding  case  we  have 

x,  '  x*  *  x3                 *■  2u,  *  2u*=  -  TZ 
and    —  i<  -t —  +-  •  •    —  _  4  : —  ■+■  : —  i  : —  -  —    H  • 

X/       XA      XS  Xh       U'    *    lri         U|    "    1V/  UX.  *    1VZ.       UJL~    1V*  9>T 

This  system  may  be  written 

fx,  t  xt  t-  x3  f  .  ....+  x„  /•  2u,  *  2u,,=  -  £ 
V—  -y  f-  —  4-  •  •  •  ™  i  —  ■+■  — - —  —  

(x<    xa.   x,  x*  r£  a>7 

and  is  now  easily  solved  for  u,  and  u^,    By  a  generalization  of  this 
principal  we  may  solve  equations  which  have  more  than  two  pairs  of  complex 
roots . 

Case  III.     Coincident  Roots. 

Coincident  roots  are  detected  in  an  equation  by  the  root- 
squaring  process  whenever  one  of  the  doubled  products  always  remains 
in  the  first  row,  or  whenever  one  of  the  coefficients  of  a  transformed 
equation  is  less  than  the  corresponding  coefficient  of  the  given  equation. 
In  such  cases  the  equation  is  broken  up  into  linear  and  quadratic 
fragments  and  solved  much  the  same  as  in  the  case  of  complex  roots, 
through  use  of  the  relationships  of  the  roots  and  coefficients  of  the 
equation. 

The  root-squaring  method  has  also  been  extended  to  the 

solution  of  transcendental  equations,  by  expanding  the  equation  into  a 
*  For  the  ratio  of  two  coincident  roots  never  becomes  <.  1. 
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Taylor  series,  neglecting  the  remainder  term,  and  then  treating  the 
resulting  polynomial  as  an  algebraic  equation, 

(f)  Bernoulli's  Method.  * 

In  1728  Danial  Bernoulli  devised  a  method  wholly  different 
in  principal  from  any  which  were  then  known.     He  made  use  of  a 
difference  equation  by  means  of  which  a  series  of  values  could  he 
calculated,  the  ratio  of  any  two  of  which  in  general  approached  as  a 
limit  the  absolutely  greatest  root  of  the  given  equation. 

(g)  Solution  by  Inverse  Interpolation. 

Newton's  formula  for  interpolation  is  as  follows: 

n(n-l)Dz   ,  n(n-l)(n-2)D,  , 
y  =  --o  ♦  nD,   +      21"  +  — 3>    

where  ya  is  the  initial  value-  of  the  function,  n  is  the  ordinal  number 
of  the  term,  and  B,  ,       ,  D  ,  etc.,  are  the  first,  second,  and  third 
differences.     If  we  wish  to  find  the  valus  of  the  argument  corresponding 
to  a  given  value  of  the  function  we  solve  the  above  formula  for  n. 

If  only  the  first  differences  are  taken  into  account,  then 

y  -  yo 


v  —  v« 

Similarly  ns  * — r-r 


d,  f  s-cn-i;r^' 

and  n  ^  D,  ^(n-^D^f  l(B-l)(n-2)D3  * 

Suppose  we  are  to  find  the  real  root  of 
3  X  +  5= 0 

which  lies  between  x=.l  and  X=.2, 

We  first  form  a  tab  Is  of  differences  of  the  function 
2S  3   -•  4x*   f  5. 

*    The  Calculus  of  Observations  -  Whittsker  and  Robinson,  Pages  98-99. 
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x  y  D,  D3 

-  1  0 

0  5  -  8 

-  3  6 
12  -2 

-  5  6 

2-  3  4 

-  1  6 

3-  4  10 

9 

4  5 

*'.'e  wish,  of  course,  tho  value  of  x  for  which  y.=  0.    Using  the  above 

derived  formulae  for  successive  approximations  we  have 

n  =  °J^i=  .4000 

1  D,  -  5 

y  -  yg     -  2     

J=    D,   tfivljD^Xiv1)^-  2)D3~  -   •  +  (-i.        ;   +  "*.3832 

 Y  -   yA  _  -   2  _ 

ny  =    ^    t  ;/,^-l)Dzt  ii  ^-Ijin^-^"  -  C  f  (-1.2336)  f  .9972"   • 3 3194 

Hence  x=x0  4  n=  1.3  8194 

Divided  differences  may  also  he  used  in  the  solution  of 

equations . 

2.    Simultaneous  Equations  in  Several  Unknowns. 

(a)     The  ITewton-Eaphson  Method  for  Simultaneous  Equations. 

The  roots  of  simultaneous  algebraic  and  transcendental 
equations  in  several  unknowns  can  be  found  by  the  Newton-Raphson  method. 
We  will  consider  the  case  of  two  equations  in  two  unknowns  since  all  the 
essential  points  are  thereby  involved. 

Let  the  given  equations  be 

/(x,y)=0,  (1) 
^     ,  =0. 

-  30  - 


r 


C 


Now  if  r.0 ,  yo  be  approximate  values  of  a  pair  of  roots  and 
h,  k  be  corrections,  so  that 

x  =  x  0  f  h, 

:/=  y«  * 

then  (l)  and  (2)  become 

/>(x„  *  /  k)=  0,  (3) 

*  b,  y,   /•  k)  =  0.  (4) 
"xpanding  (3)  and  (4)  by  Taylor's  theorem  for  a  function  of  two  variables, 
and  neglecting  higher  powers  of  h  and  k,  we  have 

rfU.  *  h,  y#  f  k)=  /(x.,  yj  *  h(|l)#   f  k(|£)o=    0.  (5) 
*  h,  y.  ,  k)=>(x#1  yj  /  h(||)o   *  k(§|).=:    0.  (6) 
Solving  these  equations  by  determinants,  we  find  the  first  corrections 
to  be 

-  y.) 


^  / 


D 


(7) 


k.  = 


•44^  yj 
'lit)    -0(x      V  ) 


(8) 


wher' 


D  = 


(dp) 


ayi 


O) 


Additional  corrections  can  be  found  by  repeated  application  of 
these  formulae  with  the  improved  values  of  x  and  y  substituted  at  each 


step.     (The  notation  (-^-^)    means  the  value  of  (5^-^    \vhen  x    and  v„ 

o 

substituted  for  x  and  y. ) 


are 


!•  3 1  — 


Focample  :    Given  the  equations 

(x  f  31oS/(>x  -  y%  0, 

(2x    -  xy  -  5x  f  1  =  0. 

By  plotting  the  graphs  of  these  equations  on  the  same  set  of 

axes,  v/e  find  that  they  intersect  at  the  points  (1.4,  -  1.5)  and 

(3.4,  2.2).  Let 

(*$(x,  y)=x  f  31og  /6x  -  y*~ 
( 

(V(x,  y)=  2x    -  xy  -  5x  f  1. 

Then  ^-=1  ,  (where  M=  0.43429) 

x  x 

-2v 

^x       4x  -  y  -  o,  jt7--  x. 
Now  since  x0=3.4,  y,=i  2.2,  we  have 

t(**,  yo)=-  0.1545,^  (xe,  y0)=  -  0.72, 
1.383,  -4.4, 
(^=6.4,  ^  -  3.4 
Substituting  these  values  in  (7),  (8),  and  (9)  we  find 

h,  =  0.157,  k,  =  0.085 
Hence  x,  — S.4  -  0.157=  3.557,  y/=  2.285 

Now  substituting  x,  and  y,  for  x  and  y,  we  get 

hAs  -  0.0685,  k^ss.  -  0.0229. 
Hence  x^=  3.4885,  y^s  2.2621 

Similarly  the  third  approximations  are 

x3  =i  3.4732,  y3=  2.25154 
♦Taken  from  Numerical  Mathematical  Analysis,  J.B.  Scarborough,  PagQ  189 
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(b)    The  Method  of  Iteration  for  Simultaneous  Equations. 
Given  the  two  equations 

/(x,  y)  =  o  (i) 
f  (x,  y)=  o  (2) 

we  first  write  the  given  equations  in  the  forms 

x=F,(x,y)  (3) 
y  =  FJL(x,y).  (4) 

Then  if  (x,  ,  y0  )  "be  the  approximate  values  of  a  pair  of  roots,  improved 

values  are  found  by  the  steps  indicated  below.- 


1st 

approximation 


2nd 

approximat  ion 


(xy=  F,(xa  ,  y0  ), 

(y,=  FJX,  »  y»  )j 

(x  =  F  (x,  ,  y,  ), 

( 

etc. 


(5) 


(6) 


Hot  every  form  of  (3)  and  (4)  will  converge  and  others  will 
converge  very  slowly,  hence  care  must  be  taken  in  the  selection  of 
these  functions.     It  can  be  shown  *  that  the  iteration  process  for  two 
unknowns  converges  when,  and  only  when,  the  two  conditions 


<5?J 

ax 

+ 

ax  1 

-r 

5y 

(7) 
(8) 

hold  for  all  points  in  the  neighborhood  of  (x4  ,  y„ )  •    Further,  in  order 
for  the  convergence  to  be  rapid  enough  to  make  the  method  advisable  in 
any  given  problem  it  is  necessary  that  these  two  values  (7)  and  (8)  be 
much  less  than  unity.    Thus  the  method  of  inter  at  ion  should  not  be 
chosen  at  all  for  some  examples. 

*    Numerical  Mathematical  Analysis  -  J.  B.  Scarborough,  Pages  193-195. 
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IV.     GRAPHICAL  METHODS  OF  SOLUTION. 

Graphs  are  much  used  in  the  solution  of  equations  when  a 
rough  approximation  is  all  that  is  desired.     A  function  f(x)-0 
may  be  solved  by  plotting  the  graph  of  y=f(x),  in  which  case  the  real 
roots  of  the  equation  are  the  points  where  the  graph  corsses  the  x-axis, 
or  a  function  f(x)=0  may  be  broken  up  into  two  other  functions  f,  (x)=f  (x) 
(one  or  both  of  which  may  be  standard  curves)  and  then  plotted  separately 
on  the  same  axes,  in  which  case  the  abscissae  of  the  intersections  of  the 
two  curves  are  the  desired  roots, 
1.    Solution  by  Repeated  Plotting. 

(a)    Functions  of  the  Form  f(x)=  0. 

The  method  of  solution  by  repeated  plotting  is  as  follows: 

(1)  Find  an  approximate  value  of  the  root  either  by  a  rough  graph,  or 
otherwise, 

(2)  Plot  on  a  large  scale  a  small  part  of  the  graph  y=f(x)  in  the 
neighborhood  of  the  desired  root.  An  additional  figure  of  the  root 
will  thus  be  found, 

(3)  Again,  plot  on  a  still  larger  scale  a  small  part  of  the  graph  near 
the  new  approximation,  getting  a  still  more  accurate  value  for  the  root. 

(4)  Continue  this  process  until  the  desired  accuracy  is  obtained. 

This  method  may  be  illustrated  by  the  following  example.* 
Find  the  positive  real  root  of 

x  -  ^,(0.7854  -  *Tfl  -  **)  _  0. 

*    Taken  from 

Numerical  Mathematical  Analysis  -  J,  B.  Scarborough,  Page  177 
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We  first  compute  the  value  of  the  left  member  for  several 
values  of  x,  as  given  in  table  (l).     This  table  shows  that  a  root  lies 
between  0.5  and  0.6.     Hence  we  plot  the  graph  of  the  given  equation 
from  x=0.5  to  x=0.6  and  assume  it  to  be  a  straight  line  within  this 
interval.     The  result  is  given  in  figure  (a)  and  it  shows  at  a  glance 
that  the  root  is  about  0.56  or  0.57.    We  therefore  compute  table  (2) 
and  plot  the  results  as  shown  in  figure  (B).     This  graph  shows  that 
the  root  is  about  0.579.     Continuing  the  process  in  this  manner  by 
computing  tables  (3)  and  (4)  and  plotting  the  results  on  still 
larger  scales  as  shown  in  figures  (C)  and  (D),  we  find  the  desired  root 
to  be  x  =0.57936  to  five  places. 


(1) 


x 
0.4 
0.5 
0.6 


f(x) 

-  0.42 

-  0.26 
*  0.14 


(3)  0.579 
0.580 


f(x) 
-  0.001 
f  0.003 


(2) 


-0.26 


0.56 
0.57 
0.53 


-  0.092      (4)       0.5793  -  0.0005 

-  0.030  0.5794  4  0.0003 
+  0.003 
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(b)    Functions  in  the  Form  f  (x)=.f  (x) 

Another  convenient  method  of  graphical  solution  is  to  write 
the  equation  f  (x)  =  0  in  the  form  f,  (x)=f^(x),  which  may  us\ially  be 
done  in  several  ways.    For  example,  the  equation  x3  -  4x  f  6=  0  might  be 
broken  up  into  y,  =•  x3,  and  y^  —  4x  -  6. 

Then  the  graphs  of  y-=f,(x)  and  y  — f^(x)  are  drawn,  the 
abscissae  of  the  points  of  intersection  of  these  two  curves  being  the 
required  roots  for  they  evidently  satisfy  the  equation  f(x)=-0.  Now 
the  steps  outlined  in  the  above  solution  are  followed  except  that  the 
portion  of  the  graph  to  be  plotted  is  in  the  neighborhood  of  the  points 
of  intersection  of  the  two  curves  instead  of  the  points  where  the  curves 
cross  the  x-axis. 
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2.     Solution  by  Lill's  method.* 

In  order  to  determine  the  roots  of  the  equation 
f(x)=anxn   f  an_,x        f  a,  x  +  a6=0 

the  coefficients  a  ,  a^  ,  afl  are  represented  by  lengths,  expressed 

in  terms  of  a  suitable  unit  l^mm.     These  distances  are  set  off 
successively,  so  that  the  initial  point  of  the  line  corresponding  to 
any  one  coefficient  coincides  with  the  terminal  point  of  that 
corresponding  the  previous  coefficient;  each  line,  moreover,  makes  an 
angle  with  its  predecessor  of  90*  measured  in  the  direction  of 
rotation  of  the  hands  of  a  clock.      The  accompanying  diagram  illustrates 

this  procedure  for  a  fourth  degree 
function,  where  for  the  sake  of  simplicity 
all  the  coefficients  are  taken  to  be 


f—U-t 

i   '> 


B, 


'f— if  ~t\ 

1   A*  1, 

B\      «.  A, 


i  >  a. 

\ 

\ 


\ 


\ 


X 


6, 


Positive.  A  A  ,  A  A  ,  etc.  represent 
respectively  the  coefficients  ay,  a^  , 
etc.    By  this  means  a  broken  line 


A  A  A  A  A  A  is  derived.     The  values  of 

9       T       7  © 


x  are  represented  by  taking  a  point  E 
8,  f  on  A/  A^  at  a  distance  1^  (the  unit  of 

length  x-wise)  and  erecting  a 
perpendicular,  and  measuring  on  this  line  the  values  of  x  for  which  the 
function  is  to  be  evaluated.    "When  x  is  positive  the  measurement  is 
made  from  E  upwards  and  vice  versa.    For  x=  x,  we  obtain  a  point  Xf 
such  that  EE.  =:   x  •  1  mm. 

*    Practical  Mathematical  Analysis  -  H.    Von  Sanden,  Pages  57  -  60. 
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Join  X/  to  A    j  then,  angle  X/A/E=y[j  ^s  such  that 

tan  f>(=  x,  . 

Let  A/X/cut  A^A    in  a  point  B/  ,  then 

-  s_  x  .  1  mm. 

Further,  B.A.=  (ax,   +  a,)-lmm. 

Drawing  a  perpendicular  at  B/  to  A/B/  we  determine  BA  its  point  of 
intersection  with  A  A  ;  then  angle  B,B/  A  =■  £  ,  and  the  distance 

A3BZ  =  *  SX<  )  *  1<s.Tmn« 

Hence  Af  BA  ~    ^%xf  *  S  x/  +  a*  )  *  1<Lmm« 

Now  B^B   X  B,  3^  a-nd  By  B^ X  \B3  ,  and  it  is  easy  to  show  that  the  distance 
By  A^  s  (avx/  a  ^x3  +  iizx*  *  a,x,  *  a,  )  •  l^mm;  that  is  to 
say  this  distance  represents  the  value  of  the  function  for  x=x;,  measured 
in  terms  of  the  unit  of  length  1  , 

To  determine  a  root  of  the  equation  f(x)  — 0,  the  value  of  x, 
and  consequently  the  position  of  E  on  the  line  Ay  A  ,  has  to  be  altered 
so  that  the  terminal  point  B^  of  the  sequence  A/B/B^BjB^  coincides  with 
A.  •     This  would  then  give 

B y Afe  —  0,  i.e.  f(x)=0. 
For  this  purpose  the  use  of  transparent  millimetre  paper  is  of  great 
advantage.     This  is  fixed  at  the  point  A,  by  means  of  a  needle,  so  that 
the  paper  may  be  rotated  about  this  point.     This  is  then  done  until  one 
of  its  lines  passes  through  Ay  and  X/  ,  and  then  the  polygon  A^B^B^B^ 
may  be  traced  without  drawing  a  single  line.     The  paper  is  now  turned  until 
a  position  is  found  for  which  B     falls  on  A  .     The  corresponding  point  x 
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is  then  marked  "by  a  prick.     It  is  advisable  to  turn  the  paper  by  a 
very  small  angle  at  a  time,  noting  the  change  in  the  position  of 

since  conclusions  regarding  the  reality  of  the  roots  may  be 
drawn  from  it.    When  the  coefficients  have  different  signs  the  follow- 
ing rule  is  to  be  observed:    draw  a  rectangle  associating  each  side 
with  an  arrow  marking  the  direction.    Beginning  with  the  coefficient 
a„  of  the  highest  power,  we  let  it  correspond  to  the  upper  horizontal 
side  of  the  rectangle,  the  coefficient  to  the  neighboring  side, 

etc.,  going  in  clockwise  direction.    Negative  values  are  set  off  in  the 
direction  opposite  to  the  arrow.     If  a  coefficient  is  zero  then  the 
lines  corresponding  to  the  preceding  and  succeeding  coefficients  lie 
in  one  line.     The  accompanying  diagram  represents  the  solution  for 


B, 


the  case  of  the  equation 

4x3  *  6xK  -  3x  -  1.9=  0. 


We  may  now  divide  the  equation  by 


x  -  x/  and  solve  the  resulting  quadratic. 


By  following    the  sequence  of  lines 


x,   of  the  given  equation  is  found  to 


have  the  value 


giving  the  polygon  A,B,B^Ar,  a  root 


X,  -  tan      =  0,67. 


By  proportional  parts  we  find  that  the 


sides  of  the  polygon  A,B,B^B 
equation  of  the  second  decree 


may  immediately  serve  for  representing  the 


•  'JL  Z  +  s.'  x  /  ;/  s  0 

■which  provides  the  remaining  roots.     These  can  be  found  by  drawing  a 
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semicircle  on  A, A*,  as  diameter.     Its  points  of  intersection  Uy  and  U 

with  B/  B^are  each  joined  to  A/  and  A^..     The  two  polygons  A/U/  A-  and 

A  TJ  A  will  now  furnish  us  with  the  remaining  roots  x.  and  x 
/  ^  j  a.  , 

and  ^    are  the  angles  between  the  line  representing  the  coefficient  a^, 
that  is  A,B,  ,  and  the  lines  A/TJ    and  A,!^  respectively.  Accordingly 

x,  =  tanrf  =    -  0.40:  x  ,=  tan/  =-  1.77 
are  the  values  of  the  roots.    Had  the  semicircle  on  A  A-.  not  intersected 
the  line  B,  B^  ,  this  would  have  been  an  indication  that  the  roots  are 
not  real, 

V.    NOMOGRAPHIC  METHODS  OF  SOLUTION 

Graphical  methods  of  solution  have  in  many  instances  been 
superseded  by  nomographic  methods,  in  which  a  diagram  is  prepared  to 
be  used  for  a  wide  number  of  cases.     The  field  of  nomography  is  not 
primarily  interested  in  the  solution  of  equations,  hence  the  following 
car  hardly  be  considered  as  typical  nomograms.    However,  they  are 
valuable  in  solving  equations,  for  they  may  be  used  over  and  over 
again  with  different  numerical  data, 

1.     Analytic  Chart  for  Solution  of  Algebraic  Equations.* 

Ye  may  use  the  rulings  of  a  sheet  of  rectangular  coordinate 
oaper  to*  solve  graphically  algebraic  equations  of  the  4th,  3rd,  and  2nd 
degrees . 

Case  I.     If  we  draw  the  parabola  y  —  2x  and  the  circle  (x  -  h)  f 
(y  -  k)  =  r    with  center  at  (h,  k)  and  radius  r,  the  ordinates  of 
their  points  of  intersection  are  found  algebraically  be  eliminating 
x  between  these  two  equations  and  solving  the  resulting  equation  for  y. 
*    Graphical  and  "Mechanical  Computation  -  Joseph  Lipka,  .  Pages  26-28. 

-  40  - 


A 


f 


From  the  first  equation  we  have  and  substituting  this  in  the  second 

equation  we  get 

if-  h  f  f  (y  -  k)\=  rL,  0r 
y¥  +  4(1  -  h)y*  -  8ky  f  4(hx  f  k*  -  r*)  »  0. 
If  we  divide  this  equation  by  t  ,  --here  t  is  an  arbitrary  number,  we  get 
(i/  +         *  hVf      k.y.       ,  (h*  fk^-  r")  _ 


or  z*  f  az2"  *  bz  f  c-  0, 

where  z  =  1  ,  a=  ilL^Jil     b  =  _  8k  4(h*  »  k*  -  Q  . 

t  t*-  t5  if 

Conversely,  the  real  roots  of  the  equation  z    f  az    *  bz  f       0  are  found 


by  measuring  the  ordinates  of  the  points  of  intersection  of  the  parabola 

y  =  2x  and  the  circle  with  center  at 

H-  4  -         k=  _  ^     and  radius  r  =\lh*  +  k*  -  °£ 
4  » 

where  t  is  an  arbitrary  number,  and  dividing  these  ordii;fetes  by  t.  The 
use  of  t  allo'v/s  us  to  throw  the  center  of  the  circle  to  a  convenient 
point  (always  to  the  right  of  0Y,  or  as  near  to  or  as  far  from  the  vertex 
of  the  parabola  as  is  convenient)  such  that  the  circle  will  not  cut  the 
parabola  at  an  angle  too  acute  for  accurate  reading  of  the  ordinates. 
Thus  the  one  parabola,  y-s  2x,  will  serve  for  finding  the  real  roots  of 
all  equations  of  the  fourth  degree.     To  solve  the  complete  equation 
v^  +  pv3  fqv    +  rv  +  s=0,  we  first  substitute  v=.  z  -  H    and  this 
equation  takes  the  form  z    +  ?zz  +  bz  t  c=0. 

"xample  1.*  Let  us  find  the  real  roots  of  the  equation  z    +  z  -  1=0. 

Here  a=  0,  b  =  1,  c=-  1.    Hence  p-=l  *nd  k  =  — -it3.     If  we  choose  t—  2, 

*    Taken  from  Graphical  and  7:echariical  Computation  -  Joseph  Lipka,  Page  27 
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the  center  is  tie  point  (l,  -  1)  and  the  radius  is    V~6  =  2.45.  The 
circle  cuts  the  parabola  in  tv/o  points  whose  oi'diriates  are  approximately 
y  =  1.4  and  y  =  -  2.4.     "ence  z  —  *■  =  0.7  and  -  I.E. 


y 


-5 


Case  II.     If  c=0  ir  the  equation  z*  f  ft»x  f  bz  +  c=  0,  this  equation 
"becomes        f  az*"  f  bz  =  0,  or  z(z3  f  az  +  b)=  0,  or  z-  0  and 
z3  +  az  f  b  —  0,     One  of  the  roots  being  zero,  the  circle  will  pass 
through  the  origin  or  vertex  of  the  parabola,  and  the  other  points  of 
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intersection  will  give  the  real  roots  of  the  cubic  z3  f  az  ■+■  b=  0. 
Hence,  the  real  roots  of  the  equation  z3  f  az  f  b=0  are  found  "by  measur- 
ing the  ordinates  of  the  points  of  intersection  of  the  parabola  y*=  2x 

4  —  at**  bt^ 

and  the  circle  with  center  at  h= — j — — ,  k-:  —  — —    (where  t  is  an 

%  8 

arbitrary  number)  and  passing  through  the  vertex  of  the  parabola,  and 
dividing  these  ordinates  by  t.     To  solve  the  complete  equation 
v    +  pv    ¥  qv  f  r  =-0,  we  substitute  v=- z  —       thus  getting  the  reduced 
cubic,  and  proceed  as  above. 

Lxample  2.*  Let  us  find  the  real  roots  of  the  equation  v3  -  3v*  \  v  +  4=  0, 

Let  z  — v  f  lj  then  the  equation  becomes  z3  -  2z  -  5—0,     Here  the  center 

i  5  3 

of  the  circle  is  at  h  =  1  f  jgt  =  3 ,  k=  -gt  =  5,  if  t  =  2.    T/'e  read  y=4.18j 
hence  z  =  J=  2.09,  and  v=  z  -  1  =  1.09. 

Case  III,     If  we  drav:  the  parabola  y*=  2x  and  the  straight  line  y=;mx  t  k 
of  slope  m  and  y-  intercept  k,  the  ordinates  of  their  points  of  intersection 
are  found  from  the  equation  y2"  —  ^y  f  ±il-o.     If  we  divide  this  by  fc* 
(where  t  is  an  arbitrary  number),  we  get 

W  ■  r3fi&  +  mT^0     or      z    +  az  *  b=  °» 

7/here  z  =  ^L,  ~,  b=.i^L.     Conversely  the  real  roots  of  the  equation 

t  mt  mt*" 

z*~  f  az  +  b=0  are  found  by  measuring  the  ordinates  of  the  points  of 

intersection  of  the  parabola  y  —  2x  and  the  straight  line  of  slope 

2  bt 

rn-  —    and  y-  intercept  k^r  and  dividing  these  ordinates  by  t. 

3."u  a 

1'ote  that  the  one  parabola  y*=  2x  will  serve  for  finding  the  real  roots 
of  all  fourth,  third  and  second  degree  equations. 

Scrapie  3.*    Let  us  find  the  real  roots  of  the  equation  z2-  -  1.45z  -  5.6=0. 

2  2 

Here  the  slope  of  the  line  is  m-  '         -  -         and  its  y-  intercept  is 

1.45t      1.45  J 

k  ~-  —  -  3.86,  if  t  -1.    We  read  z=  y--  1.75  and  3.20. 

*  Taken  from  Graphical  and  Mechanical  Computation  -  Joseph  Lipka,  Page  28 


2.    Rectangular  Chart  for  the  Solution  of  Cubic  Equations.* 

The  roots  of  the  cubic  equation  z    f  pz  f  q=0  depend  upon 

the  values  of  the  coefficients  p  and  qj     thus  the  roots  are  functions  of 

p  and  q.     If  we  choose  x  =2DCp,  y»mq  for  our  perpendicular  scales,  the 

equations  of  our  representing  z-  curves  have  the  form  7/   +  zx  f  mz  ~  0, 

a  set  of  straight  lines.     In  the  chart  (below),  these  straight  lines 

are  constructed  for  values  assigned  to  z,  viz.,  znO,     -  0.1,  -  0>2, 

 ,  -1.3,  and  lying  within  a  square  bounded  by  p~-  1  and  qn"*  1. 

Eaoh  line  is  constructed  by  means  of  two  points  on  it.     Thus  for  the 

straight  line  marked  £  =  0.3,  we  have  0.027  +  0.3p  f  q=  0,  and  to 

construct  this  line  we  may  use  the  points  for  which  p-ssl,  q—  -  0.327 

and  ?=  -  1,  q=  0.273. 

On  this  chart  we  may  read  the  approximate  real  roots  of  any 

cubic  equation  for  which  p  and  q  lie  within  the  limits  -  1  and  4-  1. 

3 

Thus  for  the  equation  z    +  0.2z  -  0.8=0  we  have  p=0.2  and  q=  -  0.3, 
and  we  read  z—0.86,  interpolating  this  value  of  z  between  the  lines 
marked  z  =  0.8  and  z  =  0.9.     According  as  the  point  (p,q)  falls  outside 
of,  on  the  boundary  of,  or  within  trie  triangular  shaped  region  on  the 
left,  we  can  read  one,  two  or  three  values  of  z,  and  the  corresponding 
•cubic  equation  has  one  root  only,  three  real  roots  two  of  which  are 
equal,  or  three  distinct  real  roots.     Thus  for  the  equation  z^  -  0.7z 
•f  0.1-0  we  have  ps-  0.7  and  q ~  t  0.1  and  we  read  z  ~  -  0.91 
f  0.76,  f  0.15. 

If  the  values  of  p  and  q  lie  beyond  the  limits  -  1  and 
f  1,  the  chart  may  still  be  used. 

*    Graphical  and  Mechanical  Computation  -  Joseph  Lipka,  Pages  35  -  37 
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Let  z=kz',  s>nd  the  equation  z3  +  pz  f  q  =  0  becomes  k3  z'    +  pkz'  f  q=0, 
i3  P. 

or  z'     +  £i2    +  \=  0.    We  may  now  choose  k  so  that  p*  and  q'   lie  within 
k«-  k3 

the  limits  -  1  and  +  1,  and  read  the  corresponding  values  of  z'  from 
the  chart.     The  roots  of  the  original  equation  are  then  z=  kz' .  Thus  to 
solve  the  equation  z*  +  z  -  2=  0,  let  z=2z'.    We  get  z'3    +  0.25z'  -0.25=0, 
for  which  we  read  z  —  0.50,  and  hence  z  =  l.CG. 


-  45  - 


c 
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In  a  similar  way  we  may  "build  a  rectangular  chart  re- 
solution of  the  quadratic  equation  b    ♦  pz  f  q  =  0,  or  for  any 
trinomial  equation  z*1  +  pz"  +  q=0. 
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SUMMAET 

In  this  thesis  four  types  of  solutions  of .  algebraic  and 
transcendental  equations  have  "been  considered,    "e  first  dealt  v/ith 
the  literal  methods  of  solution  which  cannot  in  general  be  applied  to 
equations  of  higher  degree  than  the  fourth. 

Tartaglia's  solution  of  the  cubic  and  the  trigonometric 
solution  of  the  irreducible  case  are  famili&.r  methods  of  procedure. 
Ferrari's  solution  of  the  quartic  is  also  well-known  hut  the  series 
formula  for  the  root  is  of  quit?  recent  development,  and  different  in 
character  from  the  other  literal  methods.     It  consists  of  a  series  in 
which  the  coefficients,  of  the  equation  are  substituted,  giving  an 
approximation  of  the  desired  root. 

Much  of  the  space  has  been  devoted  to  the  numerical  methods 
of  solution  since  they  are,  on  the  whole,  the  most  useful.  The 
Ruff ini-Horner  method,  the  method  of  iteration,  the  method  of  false 
position,  and  the  ITewton-Raphson  method  have  much  in  common,  and  to 
bring  out  this  similarity,  the  equation 

y*  -  2x  -  5-  G 
has  been  solved  by  all  four  methods. 

The  root-squaring  method  of  Graeffe  is  unique  in  that  it 
can  be  used  to  find  all  the  roots  of  an  equation  whether  real  and 
unequal,  real  and  equal,  or  complex.     The  given  equation  is  transformed 
into  one  v/hose  roots  are  high  powers  of  the  roots  of  the  original  equation. 
These  are  widely  separated  and  consequently  easily  found.     The  formulae 


for  inverse  interpolation  provide  an  interesting  solution  which  is 
essentially  iterative,  and  which  quickly  yields  a  high  degree  of 
accuracy.    Simultaneous  equations  in  several  unknowns  may  be  solved 
"by  the  Newton-Eaphson  method  and  the  method  of  iteration. 

Graphical  methods  often  serve  for  the  purpose  of  finding 
first  approximations  of  the  values  of  roots  hut  can  also  he  employed 
in  obtaining  close  approximations  by  the  method  of  repeated  plotting. 
Lill's  method  is  a  particularly  interesting  diagramatic  method  of 
solving  for  the  roots  of  an  equation  where  the  coefficients  are 
represented  by  lengths  and  the  root  is  read  as  the  tangent  of  an  angle 
formed  thereby. 

Nomographic  methods  are  especially  valuable  in  that  a 
construction  is  made  once  for  all  and  can  then  be  used  for  any  number 
of  cases.     Two  such  nomograms  have  been  here  presented  for  examples, 
but  this  discussion  has  by  no  means  exhausted  the  field,* 


*  Cf,  the  alignment  chart  for  the  solution  of  quadratic  and  cubic 
equations  on  page  112  of  Graphical  and  Mechanical  Computation  - 
Joseph  Lipka. 
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